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ABSTRACT 


The  phenomenon  of  a colored  sun  or  moon  can  be  explained  by  Mie’s  theory  of  large  particle  scattering, 
if  the  radius  of  such  particles  is  of  the  order  of  the  wavelength  of  light. 

Smoke  generated  by  large  forest  fires  in  Alberta,  Canada,  traveled  across  the  Atlantic.  In  many  places 
in  the  eastern  United  States  and  especially  over  large  parts  of  Europe,  the  sun,  moon  and  stars  were  seen 
blue  in  color  in  September  1950.  Using  Mie’s  theory  and  some  extinction  measurements,  the  radius  of  the 
particles  is  calculated  to  be  between  0.5  and  0.8  m and  the  concentration  175/cm3  and  127/cm3.  The  total 
amount  in  the  smoke  layer  is  about  4.7  to  6.5  X 107/cm2  column.  Three  size  spectra  are  assumed,  case  1 
deals  with  a uniform  size,  cases  2 and  3 with  a Gaussian  distribution  of  the  size  spectrum.  The  refractive 
index  for  water  droplets  is  chosen,  but  it  is  also  studied  carefully  for  other  substances.  The  radius  becomes 
larger  for  ice  at  low  temperatures  and  smaller  for  pure  smoke;  the  total  number  varies  inversely  to  the  radius, 
but  the  extinction  coefficient  remains  unchanged.  Scattering  by  spherical  particles  is  assumed,  because 
the  shape  plays  a very  minor  role.  For  anisotropic  particles,  the  extinction  coefficient  increases  by  less 
than  10  percent. 

The  intensity  distribution  of  the  solar  spectrum  is  calculated  between  3500  A and  7000  A after  passing 
through  the  smoke  layer.  The  computations  show  that  a pronounced  maximum  of  the  solar  spectrum 
occurs  around  4100  A and  4600  A in  all  the  cases  considered.  The  physiological  impression  of  such  a 
spectrum  is  bluish.  It  could  be  shown  that  the  "blue”  sun  observed  during  dust  storms  can  be  explained  in 
the  same  way  if  the  radius  of  the  quartz  particles  is  about  0.3  m and  a total  number  > 108/cm2.  Finally  the 
hypothetical  case  for  a "green”  sun  is  calculated.  The  computations  show  that  this  phenomenon  can  occur 
only  for  a narrow  Gaussian  distribution  of  particles  and  this  explains  the  rareness  of  the  phenomenon. 
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ON  THE  PHENOMENON  OF  THE  COLORED  SUN,  ESPECIALLY  THE 

“BLUE"  SUN  OF  SEPTEMBER  1950  * 

1.  COLORED  SUN  AND  LARGE  PARTICLE  SCATTERING 

The  phenomenon  of  a colored  sun  or  moon  is  not  so  unusual  as  it  seems  at  first.  Numerous  descriptions 
of  it  exist  in  meteorological  literature  and  a historical  review  is  given  by  Rodewald  (1952).  It  occurs  only 
under  certain  favorable  conditions  in  connection  with  (a)  dust  storms  near  loess*  areas  or  deserts, 
(b)  volcanic  ash  clouds  after  volcanic  eruptions  and  (c)  large  forest  fires  or  prairie  fires.  It  can  easily  be 
observed  under  artificial  conditions  when  homogeneous  water  droplets  or  aerosols  are  produced  in  large 
quantity.  The  "blue”  sun  and  moon  of  September  1950  was  observed  over  so  large  an  area  for  several  days 
that  many  reports  have  been  published. 

Before  dealing  with  the  phenomenon  of  September  1950,  it  seems  necessary  to  state  in  general  terms 
under  what  conditions  a blue  or  green  sun  or  moon  can  be  observed.  It  is  certainly  due  to  scattering  by 
large  particles,  where  large  means  that  the  radius  r of  the  particles  is  comparable  to  the  wavelength  X of 
light.  As  long  as  r ^ 1/10  X,  Rayleigh  scattering  applies,  which  states  that  the  extinction  is  proportional 
to  X“4  leading  always  to  large  values  of  extinction  in  the  blue  and  small  ones  in  the  red.  In  order  to  get  a 
blue  or  green  sun,  the  extinction  must  be  small  in  the  blue  or  green,  just  inverse  to  the  extinction  due  to 
Rayleigh  scattering.  This  occurs  only  if  r > 1/10  X,  especially  if  r « X and  can  be  explained  by  Mie’s  theory. 

Let  the  cross  section  of  extinction  be  E times  the  geometric  cross  section  of  the  particles;  E is  called 
the  efficiency  factor.  Then  the  extinction  coefficient 

7 = Tr*NE(a,  m),  (1) 

where  all  particles  are  supposed  to  have  the  same  radius  r,  and  N is  their  number  per  unit  volume.  The 
formulas  of  Mie  give  E as  a function  of  the  refractive  index  m and  the  size  parameter  a defined  by 

a = 2xr/X  = 2irrn  (2) 

with  n = the  wave  number. 

A simplified  diagram  of  the  function  E(a)  is  reproduced  in  Fig.  1,  showing  various  maxima.  Whenever 
the  extinction  curve  is  rising,  as  in  parts  a,  c and  e,  light  of  longer  wavelength  will  be  scattered  less  than 
that  of  shorter  wavelength,  as  in  Rayleigh  scattering.  But  whenever  the  extinction  is  decreasing  with 
increasing  o,  as  in  parts  b and  d,  the  opposite  will  be  true.  In  order  to  get  a blue  sun,  the  extinction  must 
be  larger  in  the  red  than  in  the  blue  and  that  means  the  particle  size  must  be  such  that  part  b is  the  effective 
region  for  the  wavelength  of  visible  light.  If,  however,  the  sun  appears  green,  the  extinction  must  possess 
a minimum  in  the  green  and  be  large  toward  the  red  as  well  as  toward  the  blue.  1 n this  case  the  radius  must 
be  such  that  the  minimum  between  b and  c is  effective. 

It  will  be  shown  later  in  Section  4.12  that  a size  spectrum  will  diminish  considerably  the  amplitude  of 
the  oscillations  after  the  first  peak.  Thus,  the  extinction  by  very  large  particles  responsible  for  parts  c, 

• Manuscript  received  for  publication  7 August  1952. 

tLoeaa  ia  a fine-grained,  wind-blown  deposit  of  alluvial  age.  It  occurs  in  the  dust-bowl  areas. 
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Fig,  L Schematic  representa- 
tion of  the  efficiency  factor  E 
for  extinction  as  a function  of 
the  size  parameter  a for  the 
refractive  index  m = L33 
(after  Holl  and  van  de  Hulst). 


d and  e will  be  more  or  less  neutral.  In  all  practical  applications,  therefore,  only  parts  a and  b should  be 
considered  and  the  particle  size  in  question  is  unique.  The  next  condition  to  be  fulfilled  is  that  the  number 
of  particles  for  which  r ~ X must  be  very  large.  This  can  be  seen  from  the  following  relation.  The  intensity 
of  solar  radiation  I*  on  a day  with  anomalous  extinction  received  at  the  grounu  is  given  by 

Ia\ r = /ox  exp  [ - Ox  + 7x«)  sec  f],  (3) 

where  /3x  is  the  extinction  coefficient  for  wavelength  X by  a normal  atmosphere  and  unit  air  mass;  7x  the 
extinction  coefficient  for  wavelength  X by  large  particles  containing  N parlicles/cm3  for  a path  length  of 
1 cm;  $ the  vertical  thickness  of  the  layer  (in  cm)  and  f is  the  zenith  distance  of  the  sun. 

The  absorbed  intensity  A is  then  given  by 

A = 1 - I A 

and 

In  A = -[In  /o  — (fi  + ys)  sec  £], 
or  putting  the  constant  values  together  in  c\  and  C2 

In  A = [&  + c\NsE(a)]  sec  f + cj.  . (4) 

From  Eq.  (4)  it  is  easily  seen  that  only  if 

ciNsE(a)  » 0, 

will  the  extinction  in  the  green  and  red  regions  be  larger  than  in  the  blue.  To  achieve  this,  the  total  number 
of  particles  Ns  has  to  be  large.  An  estimate  leads  to  values  of  107  to  108  particles/cm2  column,  if  the  proper 
values  for  /3x>  as  well  as  r = 0.5  m and  m ~ 1.33,  are  selected  for  the  region  3000  A to  8000  A (see  Table  5 
for  numerical  values  of  /3).  Thus  small  numbers  of  particles  are  unable  to  reverse  the  trend  of  atmospheric 
extinction  as  a function  of  wavelength,  and  these  figures  show  why  the  phenomenon  is  a rare  one. 

A very  instructive  graph  has  been  constructed  by  Foitzik  (1950)  for  y alone.  In  Fig.  2 the  relative 
values  of  y are  plotted  as  a function  of  X for  the  visible  spectrum,  assuming  that  7s  = 1 for  X = 5500  A. 
For  a small  radius  r = 0.1  the  extinction  decreases  with  wavelength,  the  factor  is  about  X“3.  This  cor- 
responds to  part  a in  Fig.  1.  For  r = 0.5  and  r = 0.6  n,  maximal  extinction  is  found  at  5500  A and  6000  A, 
respectively.  This  corresponds  to  the  first  maximum  (between  parts  a and  b)  in  Fig.  1.  It  will  be  seen 
later  (Section  4)  that  particles  of  this  size  cause  the  blue  sun  phenomenon.  Finally  r = 1.0  m leads  to  a 
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Fig,  2.  The  relative  extinc- 
tion coefficient  yof  water  drop- 
lets between  3000  A and  10,000 

o __ 

A for  radii  r - 0.1,  0.5,  0.6 
and  1.0  jl  on  a logarithmic 
scale.  The  values  are  plotted 
so  that  ys  = 1 at  5500  A,  the 
maximum  of  eye  sensitivity. 


WAVELENGTH  X 


distinct  minimum  around  5500  A,  and  in  this  case  the  sun  would  look  greenish.  It  corresponds  to  the  mini- 
mum between  parts  b and  c in  Fig.  1.  If  the  particle  size  increases  to  values  between  r = 1 to  1.5  ji,  the 
minimum  shifts  toward  the  red  and  the  sun  would  look  reddish.  It  must  be  mentioned  that  this  Fig.  2 
applies  to  a uniform  size,  but  the  general  conclusions  are  also  valid  for  a size  spectrum. 

From  the  foregoing  general  discussion  it  can  be  concluded  that  a colored  sun  and  moon  can  be  explained 
on  the  basis  of  Mie’s  theory  of  large  particle  scattering,  provided  that  the  total  number  of  particles  is  large, 
e.g.,  > 107/cm2  column,  and  the  specific  color  depends  on  the  radius  of  the  particles. 


2.  OBSERVATIONS  OF  TIIE  "BLUE”  SUN  DURING  SEPTEMBER  1950 
2.1.  GENERAL  OBSERVATIONS 

During  September  1950,  unusually  large  forest  fires  broke  out  in  Alberta  and  British  Columbia,  Canada. 
They  reached  their  peak  from  22  to  24  September  and  generated  a tremendous  amount  of  smoke  particles. 
The  smoke  moved  slowly  southeast  over  the  continent.  But  part  of.  it  traveled  very  rapidly  across  the 
Atlantic,  with  an  average  speed  of  about  110  to  120  km/hr.  Over  large  parts  of  Europe,  the  sun,  moon 
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and  stars  were  seen  blue  in  color  between  26  to  30  September.  Visual  observations  (as  far  east  as  Poland) 
are  reported  in  nearly  every  European  meteorological  journal.  It  is  not  necessary  to  evaluate  each  observa- 
tion here  because  this  will  be  done  by  the  U.  S.  Weather  Bureau,  but  to  mention  only  those  few  facts  that 
are  essential  for  this  study. 

The  bandwidth  of  the  smoke  layer  was  reported  as  about  1500  km  over  the  eastern  United  States. 
The  smoke  was  concentrated  over  Canada  within  a band  150  km  wide,  with  a sharp  boundary  in  the  north 
and  east  edges,  thinning  out  at  the  south  and  west  edges.  The  European  observations  indicate  a horizontal 
spread  of  about  3000  km  or  even  more.  The  layer  rose  to  greater  heights  and  grew  thicker  along  its  path, 
indicating  an  upslide  motion  as  well  as  horizontal  and  vertical  austausch.  The  smoke  in  the  eastern  United 
States  was  largely  confined  to  a layer  whose  base  was  near  8000  ft  and  whose  top  was  at  or  above  15,000  ft 
(Wexler,  1950).  Over  the  Atlantic,  pilots  reported  the  lower  boundary  at  about  17,000  ft  (5.6  km)  at 
59°  N,  25  to  30°  W.  On  26  September  over  Leuchars,  Scotland,  a pilot  observed  the  lower  boundary  to  be 
at  31,000  ft  (9.5  km)  and  the  upper  boundary  at  43,000  ft  (13.1  km).  Thus,  the  thickness  of  the  cloud 
amounted  to  about  3.7  km  and  was  situated  in  the  tropopause  region  (Bull  1951).  On  27  September  over 
Cambridge,  England,  the  lower  boundary  was  found  to  be  between  9.1  km  and  9.8  km.  Pilots  flying  in  the 
lower  part  of  the  clou 'l  reported  that  the  air  smelled  like  burnt  paper.  This  establishes  the  fact  that  the 
particles  were  smoke  particles. 

The  colors  reported  in  the  eastern  United  States  varied  very  much  from  place  to  place,  with  purple, 
blue,  violet  and  lavender  the  most  common  colors.  Observers  in  Europe,  however,  reported  the  physio- 
logical impression  of  a blue  sun  and  a blue  moon,  of  a color  comparable  to  the  quartz  lamp.  The  color  varied 
in  shade  and  intensity,  and  over  England  one  observer  saw  the  moon  changing  color  on  the  evening  of 
26  September.  At  moonrise  the  color  was  the  normal  orange -yellow;  it  later  turned  into  a sea-green  and 
finally  into  a blue.  Furthermore,  a slight  color  change  of  sun,  moon  and  stars  was  observed  from  time  to 
time  after  the  first  appearance,  where  the  blue  color  was  most  pronounced,  up  to  5 and  6 October  over  Ger- 
many. Our  conclusion  is  that  all  sorts  of  colors  have  been  seen  in  the  United  States,  but  in  Europe  only 
blue  was  observed.  This  means  that  the  particle  size  over  Europe  was  more  uniform  than  over  the  United 
States.  Color  photographs  of  the  blue  Bun  have  been  taken  at  the  Naval  Observatory,  Washington,  D.  C. 
(unpublished);  another  one  is  reproduced  by  Gebke  (1951).  No  observations  of  halo  or  corona  around  sun 
and  moon  are  reported,  and  frequently  the  absence  of  such  phenomena  is  stated  explicitly.  Such  obser- 
vations give  a hint  that  the  particle  size  must  have  been  comparable  to,  or  smaller  than,  the  wavelength  of 
white  light. 


2.2.  SPECIAL  OBSERVATIONS 

Only  three  physical  measurements  of  solar  radiation  through  the  smoke  layer  exist. 

2.21.  Solar  Intensity  Measurements  By  Dietze 

Dietze  (1951)  measured  the  total  solar  radiation,  as  well  as  the  radiation  in  the  red  region,  using 
Schottfilter  RG2.  The  red  ratio  R (R  = Ir/It * with  Ir  intensity  measured  with  red  filter  RG2,  and  It 
total  solar  intensity)  for  this  filter  is  60  percent  for  a normal  day  at  the  ground.  He  measured  74  percent 
indicating  a relatively  large  amount  of  solar  radiation  beyond  6250  A.  Dietze  concludes  correctly  from  his 
measurements  and  from  his  visual  observations  that  the  particle  size  must  have  been  of  the  order  of  the 
wavelength  of  white  light. 


Table  1.  Measurements  of  the  solar  radiation  through  the  smoke  layer  at  Gotha,  Germany 
(Dietze,  1951)  on  28  September  1951. 


Time  (true  local  time)  1137  hours 

Zenith  distance  of  sun  53°.2 

Total  intensity  of  solar  radiation  It 0.29  (ly/min) 

Turbidity  factor 12.72 

True  red  ratio  R 74% 

R — f?0;  Rq  ■*  normal  value 7% 

Intensity  of  sky  radiation  on  horizontal  surface  H 0.42  (ly/min) 

H/It 2,4 


His  data  as  given  in  Table  1 are  very  interesting.  The  total  solar  radiation  amounts  to  only  0.29 
ly/min,  the  turbidity  (in  Feussner-Duhois  units)  to  12.7,  which  is  an  extreme  value,  and  the  ratio  of  sky 
radiation  (without  sun)  H to  solar  radiation  IT  increased  to  2.4.  The  normal  value  for  this  zenith  distance 
would  be  about  0.28  or  less.  This  means  that  for  this  turbid  atmosphere  the  radiation  from  the  sky  is 
larger  than  the  direct  solar  radiation,  which  is  comparable  to  that  of  a cloudy  day. 

2.22.  Extinction  Coefficients  Measured  By  Wellmann 

Another  valuable  measurement  has  been  carried  out  by  Wellmann  (1951).  He  took  spectra  of  the  sun 
between  3800  A and  4700  A,  using  the  2-prism  spectrograph  and  the  1 m- mirror  of  the  Hamburger  Sternwarte 
at  Bergedorf,  Germany.  The  dispersion  of  the  spectrograph  is  given  as  40  A/mm  for  Hy . The  mirror  was 
dimmed  by  about  5 m (m  — stellar  magnitudes).  The  spectra  were  taken  on  27  September.  On  10  October, 
spectra  of  the  sun  were  taken  in  which  the  general  weather  situation  was  similar  to  that  of  27  September, 
but  the  smoke  layer  had  definitely  disappeared  and  the  color  of  the  sun  was  normal.  The  spectra  were  taken 
at  about  the  same  zenith  distances  of  the  sun  and  it  is  assumed  that  the  atmospheric  transmission  on  both 
days  was  similar,  except  for  the  additional  extinction  caused  by  the  smoke  particles  on  27  September.  Well- 
mann evaluated  his  measurements  to  determine  the  particle  size,  but  his  method  was  a hit  crude.  We  shall 
point  out  in  Section  4.2  wherein  our  calculations  deviate.  The  two  sets  of  data  Wellmann  took  can  be 
used  to  determine  the  attenuation  of  the  solar  radiation  due  to  the  smoke  layer. 

On  a "normal”  day,  the  intensity  of  the  solar  radiation  In  received  at  the  ground  is  given  by 

In\ r = /ox  exp  [ -ftx  sec  f],  (5) 

On  the  day  with  the  "blue”  sun,  the  intensity  of  the  solar  radiation  can  be  expressed  as  in  Eq.  (3),  i.e., 

Iaxt  = /ox  exp  [ - (0*  sec  f + y\S  sec  f )]. 

If  it  can  be  assumed  (according  to  Wellmann’s  statement)  that  the  atmospheric  transmission  was  the 
same  for  both  days,  ft  = ft  - ft  the  ratio 

7^  = exp  t — yx*  sec  f]  (6) 

iVXf 

gives  the  transmission  due  to  smoke  layer.  For  further  consideration  we  reduce  this  ratio  for  f = 0 and 
write 

/ax 


I 
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Fig.  3.  Tin*  anomalous  extinction  as  a function  of  wavelength  according  to  measurements  by  Wellmann  (Bergedorf,  Germany, 
27  September  1950). 

Curve  1:  KS  373  with  the  points  indicating  the  values  found  by  Wellmann  for  f * 67°. 

Curve  2:  KS  374  with  the  points  indicating  the  values  found  by  Wellmann  for  f * 61°. 

Curve  3:  Reduced  values  for  aenith,  f «■  0°. 

Since  Wellmann  used  differences  in  brightness  Am,  and  in  order  to  avoid  the  units  of  stellar  magnitude 
m,  we  convert  his  data  by  the  relation 

Am  = —2.5  log  (I a/ In)*  (8) 

Since  no  absolute  values  have  been  measured  we  can  write 

log  {I a/ In)  = —0.4  Am  + constant 
or 

In  {I a/ In)  * —0.921  Am  + constant.  (9) 

The  values  Am  are  taken  from  Wellmann’s  paper  and  presented  in  Fig.  3.  The  ordinate  values  on  the 
left  side  are  Am;  on  the  right  In  (/a//*)  + constant.  Thus,  it  is  seen  from  the  measurements  that  the 
observed  values  can  be  best  approximated  by  a linear  relationship  between  X and  In  {I a/ In)  for  the  spectral 
region  2800  A to  4700  A.  It  is  not  permitted  to  extrapolate  this  curve  beyond  these  limits.  Curve  3 repre- 
sents the  mean  ratio  \ti(Ia/In)  for  zenith  distance  f = 0.  The  numerical  values  are: 

For  plate  Ks  373:  ln(/A//Ar)  - constant  — 4.466  X 104  dX 

For  plate  Ks  374:  In  (Ia/In)  m constant  — 4.680  X 104  dX. 

Mean  value: 

^(Ia/In)  * constant  — 4.57  X 104dX.  (10) 
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Fig.  4.  Hie  anomalous  extinc- 
tion as  a function  of  wavelength 
according  to  measurements  b j 
Wilson  (Edinburgh,  Scotland, 
26  September  1950). 


The  difference  between  the  two  sets  cannot  be  distinguished  on  the  graph.  Thus  a mean  slope  is  used  in 
Section  4.3. 

2.23.  Extinction  Coefficients  Measured  By  JFilson 

The  most  extensive  measurement  has  been  carried  out  by  Wilson  (1951).  He  took  spectra  between 
3811  A and  6338  A,  using  the  36-in.  Cassagrin  reflector  of  the  Royal  Observatory  at  Edinburgh,  Scotland, 
with  the  2 -glass  prism  trained  on  the  attached  Hilger  universal  spectrograph.  The  dispersion  is  given  as 
132,  45  and  28  A/mm  at  Jf«,  Hp  and  Hy.  The  spectra  were  obtained  between  1520  and  1525  CMT  on  26 
September  1950,  at  a zenith  distance  of  the  sun  of  71°.  The  plates  were  exposed  between  gaps  in  a low  cloud 
deck  with  exposure  times  10,  30  and  180  seconds.  Later  on  the  presence  of  clouds  prevented  further  expo- 
sures. Normal  sun  exposures  were  obtained  on  26  October  and  15  November  1950,  with  a zenith  distance 
of  the  sun  of  69°  and  73°,  respectively.  The  formalism  of  the  preceding  section  can  also  be  applied  to  Wilson’s 
observation,  since  he  uses  stellar  magnitudes.  The  extinction  coefficient  of  the  smoke  layer  for  vertical 
incidence  is  given  by 

0.921  , 

7x  *■  : Aw*, 

i sec  f 

which  follows  from  Eqs.  (7)  and  (8).  The  results  for  y\s  are  represented  in  Fig.  4.  Between  6338  A and 
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4028  A,  Wilson  considers  the  values  consistent  and  reliable  (crosses  in  Fig.  4).  Below  4028  A,  less  weight 
should  be  given  to  the  numerical  values  (open  circles  in  Fig.  4)  because  they  cover  a dpectral  range  that  was 
never  meant  to  be  utilized  in  the  instrument  employed  (instrumental  absorption,  spectrum  no  longer  in 
sharp  focus  and  other  reasons  are  listed  by  Wilson).  A mean  curve  has  been  drawn  through  the  experi- 
mental points.  The  extinction  coefficients  show  a pronounced  minimum  around  4300  A,  which  falls  within 
the  range  considered  reliable.  The  extinction  increases  toward  the  ultraviolet  and  green  parts  of  the 
spectrum;  however,  the  increase  toward  the  ultraviolet  is  too  steep.  As  can  be  seen  from  the  figure, 
Wellmann’s  extinction  coefficients  are  different,  indicating  that  the  particle  size  was  different  for  both  obser- 
vations. 

Wilson  also  evaluated  his  data  and  his  method  will  be  explained  later  in  Section  4.2. 


3.  CHEMICAL  COMPOSITION  OF  SMOKE  PARTICLES 

Samples  from  the  smoke  layer  have  not  been  analyzed,  but  only  two  possibilities  exist  regarding  the 
composition  of  the  particles. 

(a)  They  consist  of  nearly  pure  carbon  coated  with  water  or  ice. 

(b)  They  consist  of  organic  compounds,  especially  wood  oil  (suggested  by  Elsley  (1951)). 

Let  us  consider  the  first  alternative.  Tiny  soot  particles  have  definitely  been  produced  in  large  quantities. 
Because  they  are  hygroscopic  and  because  they  traveled  several  days  in  the  troposphere,  the  surface  is  wet.* 
Over  Europe  the  temperature  of  the  layer  was  around  —50°  to  — 60°C  and  therefore  the  outer  film  consists 
of  ice.*  Thus,  it  seems  fairly  certain  that  the  soot  particles  are  covered  with  an  ice  film,  the  thickness  of 
which  is  unknown. 

The  second  alternative  is  that  the  particles  consist  of  a complex  mixture  of  oils  (wood  oils)  which  are 
always  produced  in  such  a combustion  process,  together  with  wood  creosote  and  some  carbolic  acid  or  phenols. 
The  oils  leave  the  fire  area  in  a gaseous  phase,  condense  in  the  atmosphere,  coagulate  and  later  on  — in  the 
cold  upper  troposphere  — may  be  solidified.**  Since  these  oils  are  not  solvent,  coating  with  ice  or  water  is 
not  possible.  Because  they  are  volatile,  they  may  produce  the  odor  of  burning  wood  which  has  been  noticed 
by  the  pilots.  Elsley  states  that  airplanes  flying  through  the  smoke  layer  in  the  Great  Lakes  region  "were 
found  tp  be  covered  with  a thin  film  of  a pale  brown,  oily  smelling  substance,  that  was  probably  a resinous 
distillate."  Unfortunately  this  deposit  has  not  been  analyzed  chemically.  This  fact  also  points  very 
strongly  to  the  second  alternative.  After  the  particles  become  solid,  coagulation  ceases  and  particle  size 
remains  more  or  less  constant.  This  fact  seems  to  explain  the  various  colors  observed  in  the  United  States, 
where  the  size  of  these  liquid  droplets  varied  with  time  apd  place,  and  the  more  or  less  uniform  blue  color 
observed  over  Europe,  where  the  particles  were  solidified  and  changed  very  little  in  size. 

From  this  discussion  it  seems  fairly  certain  to  me  that  the  cloud  consisted  mainly  of  solid  oil  particles. 

• In  an  atmosphere  where  the  relative  humidity  is  > 70  percent,  practically  all  particles  possess  a droplet  character.  It  seems 
certain  to  me  that  the  amount  of  water  vapor  in  the  troposphere  has  been  sufficient  to  wet  the  surface  of  the  particles, 
t The  lowest  freezing  temperature  of  supercooled  water  droplets  in  the  atmosphere  is  given  as  about  — 40°C. 

••  It  is  possible  to  supercool  such  oil  droplets,  but  the  temperatures  at  which  this  occurs  are  unknown. 
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4.  RADIUS  AND  CONCENTRATIONS  OF  SMOKE  PARTICLES 


4,1.  THEORETICAL  CALCULATIONS 

In  order  to  evaluate  Wellmann’s  and  Wilson’s  measurements,  the  variables  as  well  as  the  constant  factors 
are  carefully  selected.  Some  are  not  properly  known.  Therefore,  the  errors  due  to  such  an  uncertainty 
are  studied  in  detail.  The  factors  are: 

(a)  The  size  spectrum;  two  cases  will  be  considered  in  Sections  4.11  and  4.12. 

(b)  Because  the  smoke  particles  are  probably  wood  oil,  the  refractive  index  will  be  rn  « 1.46,  at  any 
rate  1.25  < m < 2.0.  The  numerical  calculations  will  be  carried  out  for  m = 1.33,  solely  because 
the  efficiency  factor  E is  best  known  for  thisr  value  of  m.  This  question  will  be  thoroughly  discussed 
in  Section  4.13. 

(c)  A spherical  shape  will  be  assumed,  because  E is  best  known  for  this  shape.  It  will  be  shown  in 
Section  4.14  that  the  shape  factor  does  not  alter  the  results  remarkably. 

(d)  In  calculating  TV,  the  thickness  of  the  layer  is  assumed  to  be  s = 3.7  km,  according  to  the  pilot 
observations  over  Scotland  (Bull  1951).  Moreover,  it  is  assumed  that  the  particles  are  evenly 
distributed  within  the  layer.  These  two  assumptions  enter  only  into  the  calculation  of  TV.  The 
total  amount  A r»  independent  of  this  assumption,  will  be  used  for  all  further  calculation  in  Sections 
5 and  6. 


4.11.  Uniform  Size 

As  a first  approximation,  it  is  assumed  that  all  particles  possess  a uniform  size,  with  radius  r.  They 
have  uniform  refractive  index  m and  are  spheres. 

From  Eq.  (10),  it  follows,  after  differentiating  with  respect  to  X,  that 

d [In  (Ia//v)] 


and  Eq.  (7),  it  follows  that 


d\ 

d [In  (/a//v)] 
dX 


= —0.457  X 10s 


= - s(dy/dX ). 


The  left-hand  side  of  Eqs.  (12)  and  (13)  are  equal,  thus 


1 


dy/dX  = 0.457  X 106-- 


and  with  s = 3.7  X 105  cm 


dy/dX  = 0.1235 


[=]• 


Later  on  we  like  to  use  the  wave  number  n instead  of  X thus 

dy/dn  = ~ 


(12) 

(13) 


(14a) 


(14b) 


The  extinction  coefficient  y is  defined  by  the  Mie  theory  (Eq.  (1))  in  Section  1.  Differentiating  Eq.  (1) 
with  respect  to  n, 


I 


I 


18 


WAVE  NUMBER  n 


Fig.  5.  The  efficiency  factor  E is  a function  of  the  size  parameter  a f or  m -»  4/3  (after  Holl)  in  logarithmic  scale,  for  case  1. 
The  function  E is  represented  only  between  a » 2.8  and  10  (curve  1).  The  straight  line  (curve  2)  on  the  left  represents  the 
measured  values  of  y as  a function  of  n (after  Wellmann)  for  the  region  3800  A to  4700  A,  also  in  logarithmic  scale.  (7  in 
a relative  scale,  because  only  the  gradient  dy/dn  is  actually  measured.)  The  beat  fit  of  Wellmaxw’t  data  with  the  theoretical 
curve  is  given  as  a dashed  line  (curve  3);  the  arrows  indicate  the  limits  of  Wellmann Ys  data. 

dy/dn  = 1 rrW — ~~  = 2w2r*N  dE/ dor  (15) 

art 

Taking  the  logarithms  on  both  sides 

In  (dy/dn)  * In  (dE/da)  + In  (2v2r*iV)  = In  (dE/da)  + constant.  (16) 

Because  r and  N are  regarded  as  constant  within  the  layer,  the  second  term  is  a constant*  Furthermore, 
for  a given  X the  logarithms  of  the  gradients  are  equal  except  for  this  additional  constant  which  is  independent 
of  X.  This  expression  (Eq.  (16))  allows  the  calculation  of  r in  a very  simple  way.  If  dy/dn  is  known 
experimentally  (in  our  case  for  the  region  3800  A to  4700  A by  Eq.  (10)),  and  if  dE/da  is  known  from  theo- 
retical calculations  for  the  specific  value  of  m,  then  these  two  gradients  can  be  compared  by  a simple  graphi- 
cal method  suggested  to  me  by  E.  Inn.  On  one  logarithmic  graph  paper  y is  plotted  as  a function  of  m; 
on  another,  E ia  plotted  as  a function  of  a.  Thus  the  logarithm!  of  the  gradients  dy/dn  and  dE/ da  are  known 
over  a certain  range  of  n and  a.  Both  graph  papers  are  now  superimposed  and  the  beat  agreement  of  both 
curves  found,  then  this  position  of  the  two  papers  ia  fixed.  A certain  value  of  n corresponds  to  a certain 
value  of  a,  this  can  be  seen  easily  on  Fig.  5,  which  is  s composite  figure  of  the  two  graph  papers.  For  example, 
for  a - 8,  n m 2.4  X 104. 

Unfortunately  for  our  case,  Wellmann  measured  dy/dn  only  for  a small  region  of  n.  For  a first  approxi- 
mation, however,  it  is  sufficient  to  find  a values  to  corresponding  values  of  ».  The  better  approximation 
was  carried  out  numerically  by  calculating  dy/dn  end  4E/da  km  10  points  and  the  fine)  determination  of 
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corresponding  values  of  a and  n was  based  on  the  best  agreement  of  Eq.  (16).  The  radius  r is  found  from 
Eq.  (2)  by  inserting  the  values  a and  n obtained  by  the  procedure  just  described. 

Having  found  r,  the  next  step  is  to  calculate  N from  Eq.  (15) 


N = 


1 (dy\  1 
2rV \dn ) (dE/da)  * 


(17a) 


or  in  our  particular  case,  using  y for  the  region  3800  A to  4700  A as  measured  by  Wellmann  we  get,  using 
Eq.  (14b), 


0.247  da 
aV  dE 


(17b) 


Since  da/dE  is  negative  in  the  region,  N is  a positive  value. 


4.12.  Size  Spectrum 

The  last  chapter  refers  to  smoke  particles  of  uniform  size.  This  is  certainly  only  a first  approximation, 
because  the  particles  generated  by  a forest  fire  will  spread  over  a more  or  less  extended  size  range.  Since 
the  smoke  traveled  several  thousands  of  miles  from  Canada  to  Europe,  it  is  understandable  that  the  size 
spectrum  within  a vertical  column  of  the  smoke  layer  became  narrower  with  increasing  distance  from  the 
source.  The  size  spectrum  will  probably  vary  in  a horizontal  direction  along  the  smoke  layer,  as  indicated 
by  the  various  colors  observed  in  the  United  States  (Wexler  1950).  But  if  measurements  like  Wellmann’s 
are  carried  out,  one  is  confronted  with  a size  distribution  within  a column  at  a given  instant.  It  seems  very 
probable  to  represent  it  by  a Gaussian  distribution,  because  the  layer  consists  of  particles  of  the  same  chemical 
composition.* 


The  Gauss  error  function  is 


y = 


h 


(18) 


where  h is  the  precision  measure  (indicating  the  steepness  of  the  function), 
u = hx 

y = -7=  e"“2  = h<t>(u), 

V* 


where 


4>(u) 


1 

y/v 


This  can  be  normalized  by 


(19) 


If  the  distribution  of  particles  is  considered  it  can  be  written  as 


Nr  — Ch<f>(u)  = C — r=  exp  [ — h 2 log  (r/r0)2] 
V * 


(20) 


where  C is  an  arbitrary  constant,  instead  of  as  in  Eq.  (19).  The  quantity  iVr0  is  the  ordinate  at  the  peak, 


• A comparison  of  the  expected  distribution  can  be  made  with  distribution  functions  obtained  from  electron  microaoope  pic* 
tures  of  soot  particles  from  flames  (Pepperhoff  1951)  which  show  a very  distinct  Gaussian  distribution.  No  comparison  should 
be  made  with  distribution  functions  of  dust  particles,  however,  because  “atmospheric”  dust  consists  of  particles  of  different 
chemical  composition  and  different  origin.  They  are  inhomogeneous  in  every  respect  and  this  leads  to  very  irregular  distribution 
functions.  Gaussian  distribution  functions  have  also  been  obtained  for  raindrop-size  spectra,  which  are  homogeneous  particles 
in  the  sense  just  used. 
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Nr  « N,o  exp  [ - fc*  log  (r/r0)*]. 

The  total  number  of  particles  to  be  found  in  1 cm*  is  then 

rdr  C 00 

IV,--  / Nr  d(ln  r). 

T J — qo 

Putting  Eq.  (20)  into  Eq.  (21)  yields 

Nj*  - N„  f+"ex p [-  h'  log  (r/r0)*  d(ln  r)], 
J — 00 


(21) 


with 

It  follows  that 
because 


Nv  = N, 


U * h log  (r/r0). 

V*  i /*+w  j i 

A log  e 


f®  A log  e \/* 

+06 


/+« 


du 


(22) 


4>(u)du  ■=  1. 


The  efficiency  factor  EM  for  a size  spectrum  can  be  defined  and  it  has  been  defined  by  several  workers  in 
this  field.  Since  we  shall  use  the  numerical  values  of  Em  calculated  by  Foitzik,  we  shall  follow  his  derivation. 

The  extinction  coefficient  7 was  defined  by  Eq.  (1)  for  a uniform  size  r 


For  a size  spectrum  it  is 


7 1*  - *rWrE(a)  = — iVrE(a). 

4 m2 


■Vvx  = TN*  = ^fcmxwdc. 


(23) 


Let  us  now  define  a new  efficiency  factor  Em  {a)  for  a size  spectrum  following  Foitzik’s  definition  (1950), 

a* 


E*(a). 


(24) 


We  use  7^  firom  Eqs.  (23)  and  (24)  and  find  that 

J a *NrE(ct)da  - c?NuEu(a), 

or  with  Nm  from  Eq.  (21) 

fei *NrE(a)da  - ctEM(a)J'Ntd(\nr). 

Dividing  both  sides  by  Nu  yields  Nr/Nr9  **  y%  i.e.,  Eq.  (18), 

y cfyE(a)da  - a*Ejf(ot) ^ yd(lnr). 


Thus 
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Em(«) 


1 J^ofyE(a)da 

“s  fTyd(lnr) 


(25) 


Foitzik  calculated  Eu(ct)  for  a refractive  index  m = 1.33  and  for  three  values  of  h — namely,  hi  = 
7.07,  hi  = 5.196,  hi  = 3.873  — which  correspond  to  the  standard  deviation  = 0.1,  = 0.136  and  <n  = 


0.182. 


Unfortunately,  he  used  Stratton-Houghton’s  value  for  E(a ),  instead  of  Holl’s  or  Houghton-Chalker’s 
(1949)  which  are  better.  For  a first  approximation,  however,  both  E(aY s agree,  but  beyond  the  second  maxi- 
mum Stratton-Houghton’s  values  are  unreliable. 

This  definition  of  and  Em  (a)  allows  the  calculation  of  r<>  and  Nm  by  inserting  them  instead  of  y 
and  E into  Eq.  (16).  The  procedure  to  obtain  values  for  r0  and  Nm  is  analogous  to  the  procedure  described 
in  the  last  section. 


4.13.  Refractive  Indices  1.25  < m < 2.0 

Since  the  refractive  index  m = 1.33  — i.e.,  the  refractive  index  for  water  droplets  — was  used,  the 
question  arises  whether  or  not  this  is  a good  approach.  The  refraction  index  m of  the  particles  is  unknown, 
but  may  be  estimated  for  the  two  kinds  of  particles  discussed  in  Section  3. 

For  amorphous  carbon  m is  relatively  high,  i.e.,  m = 1.97;  for  water  it  is  1.33;  for  ice  it  is  nearly  the 
same,  i.e.,  m = 1.30,  and  for  ice  at  very  low  temperatures  it  drops  tom  = 1.25.  If  the  particles  consist  of 
two  spheres  of  different  refractive  indices  — e.g.,  the  inner  core  with  m = 1.97,  the  outer  sphere  with  m = 
1.25  — only  the  refractive  index  of  the  outer  sphere  has  to  be  considered,  provided  that  the  outer  sphere 
has  a certain  depth.  The  numerical  value  is  not  known  for  this  case  but  can  be  judged  as  probably  10"2 
to  10“ 4 of  the  total  radius  to  be  sufficient.*  Thus,  it  is  concluded  that  for  this  case  m will  be  close  to  that  of 
ice,  say  around  m = 1.3,  and  the  results  differ  only  slightly  from  those  obtained  for  m = 1.33,  for  which  case 
the  calculations  are  carried  out. 

Let  us  turn  now  to  the  second  alternative  mentioned  in  Section  3.  For  a mixture  of  wood  oils  the 
refractive  index  depends  on  the  substances,  but  they  all  have  practically  the  same  value.  The  range  is  about 
m = 1.45  to  1.47,  and  a mean  value  of  m = 1.46  is  perhaps  reasonable.  Thus  it  is  estimated  that  the 
chosen  value  m = 1.33  is  not  far  from  the  true  value.  Table  2 shows  the  various  values  which  might  be 
appropriate. 

The  m < 2 is  clearly  demonstrated  by  Wellmann’s  data.  Namely,  if  the  particles  are  uncoated  they 
should  absorb  strongly.  For  such  particles  the  fluctuation  of  the  E function  beyond  the  first  maximum 
becomes  less  pronounced  and  it  would  mean  a more  or  less  neutral  extinction,  which  is  not  observed.  Thus, 
this  points  very  strongly  to  the  chosen  value  of  m. 

The  value  m *=  ®,  often  ascribed  to  dust,  is  impossible  according  to  the  foregoing  reasons  for  m. 
Moreover,  it  is  excluded  because  the  extinction  would  be  constant  after  the  first  maximum.  Wellmann 

• Aden  and  Kerker  (1951)  havo  shown  that,  for  a particle  with  an  ice  core  and  a water  sphere  surrounding  it,  a thickness  of 
the  water  film  one-fifth  of  the  total  radius  or  more  is  sufficient  to  have  the  same  reflectivity  as  a pure  water  particle,  if  the 
wavelength  X — 10  cm  it  used.  The  wavelength  of  the  radiation  is  2 X 10*  higher  in  their  case  than  in  our  case.  Therefore,  a 
smaller  thickness  will  give  the  same  effect. 
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Table  2.  Value*  of  refractive  indices  (for  Sodium  Ddines). 


Ice  at  very  low  temperature 

Ice  at  04C 

Water 

Wood  oil 

Quarts 

Coal,  amorphous 


m 

1.25 

1.31 

1.33 

1.46 

1.55 

1.97 


argued  against  m = <»  only  from  the  point  of  neutral  extinction,  but  the  actual  values  of  ro,  as  discussed 
above,  are  already  the  best  argument  why  such  a value  can  never  be  considered. 

It  has  been  shown  by  van  de  Hulst  (1949)  that  for  1 < m < 2 the  curves  E against  a (m  — 1)  are  very 
similar  and  for  values  of  m closely  related,  the  E values  are  identical,  for  all  practical  purposes,  provided 
that  the  parameter  a(m  — 1)  is  used  instead  of  a.  Thus  the  E value  is  the  same  for  all  a(m  — 1)  if  to 
varies  between  1 and  2.  We  may  write 


or,  if  E is  constant, 


E[a(m  — 1)]  ~ E[ot„(m*,  — 1)], 
a(m  — 1)  = aw  — 1). 


(26) 


The  index  w denotes  the  values  for  water  droplets.  If  Eq.  (2)  is  inserted, 


(27) 


This  applies  to  the  uniform  size  as  well  as  for  a size  distribution.  That  means  that  only  one  set  of  E values 
for  a particular  m has  to  be  calculated  and  the  corresponding  values  E(a)  for  another  m can  be  obtained  by 
a simple  transformation.  This  applied  to  E as  given  by  Holl  and  to  Em  as  given  by  Foitzik.  The  value 
N changes,  also.  From  Eq.  (17a)  it  follows  with 

from  Eq.  (26),  that 

where 

* 2tV„»  \dn)  \dE/ 

But  if  Eq.  (28)  is  inserted  into  Eq.  (1)  it  follows  that 

7 * TT*NE[a(m  — 1)]  * (m,  — 1)]  (29) 

remains  unchanged  whatever  value  of  m is  used.  This  also  has  a very  important  effect  on  our  calculations. 
We  shall  get  other  values  for  r and  N (see  Table  3)  for  a given  experimental  set  of  yx,  but  all  calculations 
using  7 — i.e.,  Sections  5 and  6 — remain  unaltered  whatever  refractive  index  is  provided  (1.25  < m < 2). 
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Table  3.  Value*  of  r and  N as  functions  of  m. 

Relative  values  r/rw  and  N/Nw  are  calculated,  using  mw  as  standard. 


m 

1.25 

1.31 

1.33 

1.46 

1.55 

2 

r/r. 

1.32 

1.06 

1 

0.72 

0.60 

0.33 

N/N. 

0.57 

0.88 

1 

1.94 

2.79 

9.00 

4.14.  A nisotropic  Particles 

It  has  been  assumed  that  the  particles  are  spherical.  It  is  realized,  however,  that  this  may  not  be  the 
actual  case.  From  theoretical  calculation  Gans  concluded  that  the  anisotropy  is  of  very  minor  influence  as 
long  as  the  refraction  index  is  small,  i.e.,  m < 2 (see  Linke  1942,  p.  213).  The  values  of  E will  change 
very  little,  not  exceeding  10  percent.  Van  de  Hulst  calculated  an  increase  of  6 percent  for  cylinders  with 
m = 1.25  and  m = 1.50. 

The  increase  in  intensity  is,  of  course,  not  so  easy  to  determine,  because  the  increase  occurs  in  the 
exponent  of  Eq.  (3),  but  from  a few  numerical  calculations  the  author  is  convinced  that  the  influence  is  of 
very  minor  importance,  regarding  the  other  uncertainties  like  size  spectrum  and  m. 

4.15.  Thickness  of  Layer 

Whereas  the  thickness  of  the  smoke  layer  along  its  path  was  not  uniform,  as  can  be  concluded  from  the 
observations  on  the  ground,  the  altitude  range  and  the  order  of  magnitude,  as  found  by  the  pilot  over  Scot- 
land (Bull  1951),  will  not  have  changed  on  its  way  from  Scotland  to  Hamburg. 

The  total  number  of  smoke  particles  IV  r = Ns  is  obtained  directly  from  Wellmann’s  measurements 
and  this  value  alone  is  the  important  factor  entering  all  calculations.  The  factor  IV,  however,  is  calculated 
to  indicate  the  average  concentration  if  s = 3.7  km.  Moreover,  N is  certainly  smaller  at  the  lower  and  upper 
boundary  of  the  layer  than  within,  but  no  estimate  of  the  vertical  distribution  of  N can  be  given. 


4.2.  WELLMANN’S  AND  WILSON’S  METHODS 


Wellmann’s  determination  of  r differs  from  ours.  He  used  only  the  slope  of  the  extinction  coefficient 
dy/dn  and  determined  the  point  where  the  tangent  to  the  curve  E(a)  has  the  same  slope,  namely,  dy/dn  = 
dE/da.  Since  the  function  E(a)  is  not  linear,  his  method  will  yield  only  a first  approximation.  We  shall 
see  in  the  next  section  that  his  approximation  is  fairly  good. 

Wilson  used  another  method.  He  obviously  could  not  fit  his  data  to  an  extinction  caused  by  particles 
of  uniform  size.  Therefore,  he  assumed  that  an  overall  reddening  takes  place  in  the  atmosphere,  and  he 
represents  this  by  an  (l/X)x  law.  Thus,  he  wrote  the  extinction  coefficient 


(30) 


I 


The  first  factor  is  the  usual  extinction  by  large  particles  (see  our  formula  Eq.  (1)),  the  second  term  corrects 
for  what  he  calls  atmospheric  reddening. 

Computations  are  carried  out  by  trial  and  error  to  obtain  the  best  fit  for  all  observations  and  for  those 
between  4028  A and  6338  A.  No  fit  could  be  obtained  for  all  observations,  but  a good  fit  was  obtained  for 
this  limited  range.  The  numerical  values  obtained  are 


7 s 


0.71 

2.5  log  e 


E + 1.51  X 10" 4 [2a  (m  — l)]4  [• 


(31) 


From  2a (m  — 1),  r is  calculated  according  to  Eq.  (2)  and  the  total  number  from  the  constant  factor  in  Eq. 

(31) 

0.71  - xrWr. 


The  computed  values  agree  approximately  with  our  results  (see  Table  4).  The  assumption  of  an  overall 
reddening,  however,  cannot  be  easily  explained  and  he  does  not  attempt  to  verify  it.  This  factor  is  similar 
to  the  Rayleigh  scattering,  decreasing  rapidly  with  increasing  wavelength.  This  extinction  could  be  caused 
only  by  particles  smaller  than  1/10  X,  say  r < 0.05  They  must  be  a characteristic  part  of  the  smoke  layer 
and  be  extremely  numerous.  There  is  nothing  known  to  substantiate  his  assumption,  and  I think  it  is  much 
more  conceivable  to  explain  his  data  by  a size  distribution  of  the  particles,  instead  of  the  uniform  size 
which  he  assumes. 


4.3,  NUMERICAL  CALCULATIONS  AND  DISCUSSION  OF  RESULTS 

The  methods  described  above  in  principle  have  been  employed  to  Wellmann’s  and  Wilson’s  data.  The 
following  two  cases  will  be  considered:  case  1,  dealing  with  a uniform  size  follows  from  Section  4.11  and 
cases  2 and  3,  with  a Gaussian  distribution  of  the  size  spectrum  as  explained  in  Section  4.12. 

In  Fig.  5,  E is  plotted  as  a function  of  a in  a logarithmic  scale,  according  to  the  best  theoretical  values 
calculated  by  Holl  (1948).  In  the  same  figure  y is  plotted  as  a function  of  n in  the  spectral  region  measured 
(curve  2),  also  in  a logarithmic  scale.  In  such  a coordinate  system  dy/dn  is  then  a straight  line.  It  is 
obviously  seen  that  only  a small  part  of  the  corresponding  E(a)  can  be  compared  with  y(n). 

The  only  region  where  a good  fit  can  be  accomplished  is  between  a = 6.7  and  a = 9.6.  The  slope  of 
E(a)  beyond  the  second  maximum  at  a ^ 16  and  any  following  maximum  (see  Fig.  1)  is  too  steep  and  has 
to  be  excluded,  because  no  good  fit  can  be  found  for  Wellmann’s  data.  The  best  fit  is  indicated  by  curve  3, 
which  has  been  obtained  by  the  method  described  in  Section  4.11.  For  this  fit  dy/dn  = —0.15,  whereas 
Wellmann  observed  —0.12.  The  numerical  results  are  given  in  Table  4 (a). 

For  case  2 the  best  fit  has  been  obtained  with  Em  using  h = 7.07  or  a standard  deviation  of  <r  = 0.1. 
Higher  standard  deviations  gave  EM  functions  which  are  much  flatter  after  the  first  maximum,  the  largest 
gradient  was  —0.6  and  -0.38,  instead  of  -1.0  as  required  by  Wellmann’s  observation.  This  fit  is  much 
better  than  in  case  1,  as  can  be  seen  from  Fig.  6. 

The  numerical  results  as  published  in  Table  4 (a)  are  in  one  sense  surprising.  The  rough  comparison  of 
the  slope,  which  Wellmann  carried  out,  already  leads  to  the  correct  order  of  magnitude  for  r and  iVr.  The 
differences  between  cases  1 and  2 are  minor.  Whereas  r increases  only  by  1 percent,  N increases  by  about 
25  percent.  This  means  that  r can  be  calculated  very  exactly  by  simple  assumptions. 
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Table  4.  Radius,  total  number  and  number  density  of  smoke  particles. 

(a)  For  m - 4/3 


Place 

Hamburg 

Kdinburgh 

Author 

Wellmann 

Penn 
Case  1 

tdorf 

Case  2 

Wilson 

Penndorf 
Case  3 

Radius  r«  in  p 

0.53 

0.522 

0.52a 

0.820 

Total  number  of  smoke  particles  NT/ cm*  column 

7.5  X 10^ 

5 .2  X 10* 

6 .5  X 10* 

■ 

4.7  X 10* 

Number  density  of  smoke  particles  N (1/cm*)  . 

140 

175 

■ 

127 

(b)  For  1.10  <m  <2.0 


Hamburg 

Edinburgh 

m 

xtr^u 

Penndorf 

Wilson 

Penndorf 

w eiimann 

Case  1 

Case  2 

Case  3 

1.10 

r« 

.... 

.... 

.... 

2.60 

.... 

1.25 

ro 

Nt 

0.687 
3.0  X 10T 

0.660 

3.7X10’ 

1.04 

1.08 

2.68  X 107 

1.31 

r# 

Nt 

0.553 
4.6  X 107 

0.560 
5.7  X 107 

.... 

0.87 

4.14  X 107 

1.46 

ro 

Nt 

0.375 
1.01  X 10* 

0.380 
1.26  X 10‘ 

.... 

0.59 

9.12  X 10* 

1.50 

ro 

Nt 

.... 

.... 

0.52 

2.7  X 107 

0.54 

1.08  X 10* 

1.55 

0.313 
1.45  X 10* 

0.317 
1.81  X 10* 

.... 

0.49 

1.31  X 10» 

2.0 

1 

0.2 

6.0  X 10* 

0.172 
5.2  X 10* 

0.174 
5.8  X 10* 

0.26 

0.27 

4.23  X 10* 

Not*:  Case  1 * uniform  sise;  case  2 and  case  3 ■ Gaussian  distribution  of  sise  spectrum. 


Wilson’s  measurements  cover  a much  wider  spectral  range  than  Wellmann’s,  even  if  the  data  below 
4028  A are  excluded  (see  Section  2.23).  If  we  treat  those  data,  no  fit  can  be  accomplished  with  the  efficiency 
factors  E for  uniform  size.  However,  they  fit  the  efficiency  factors  for  a narrower  size  distribution  than 
Wellmann’s  very  well,  i.e.,  for  a size  distribution  with  <r  m 0.136  (case  3).  This  can  be  seen  from  Fig.  7, 
where  the  experimental  curve  d agrees  with  curve  b.  The  abscissas  for  both  curves  are  so  arranged  that  they 


EFFICIENCY  FACTOR 
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WAVE  NUMBER  n 


Fig.  6.  The  efficiency  factor 
Ey  as  a function  of  the  size 
parameter  for  at  = 4/3  in  a 
logarithmic  scale(after  Foitxik) 
for  a particle  size  distribution 
for  case  2.  The  function  E is 
represented  only  for  a between 
2 and  10.  The  best  fit  of 
Wellmann's  data  is  shown  as 
a dashed  line;  the  arrows  in- 
dicate the  limits  of  his  data. 


correspond  to  the  best  fit.  Foitzik’s  curve  for  a > 12  is  probably  too  flat;  there  should  be  a slight  secondary 
maximum  around  a « 13  or  14. 

This  solution  seems  to  me  more  satisfactory  than  Wilson’s  assumption  of  a uniform  size  plus  a general 
extinction  according  to  a X'4  law,  which  has  never  been  found  otherwise  in  such  studies.  It  seems-  much 
more  realistic  to  assume  a size  spectrum  to  explain  his  data. 

The  results  of  the  calculation  are  contained  in  Table  4.  The  radius  r0  = 0.820  p,  a value  which  is  a little 
higher  than  Wilson’s.  The  total  number  of  particles  Nt  is  4.7  X 107,  which  is  an  order  cf  magnitude  higher 
than  Wilson’s,  because  he  assumed  the  additional  small  particles.  A comparison  of  both  observations  reveals 
that  some  significant  discrepancies  exist.  Although  the  number  of  particles  — comparing  cases  2 and  3 — 
are  about  equal,  the  particle  sixes  are  not. 

The  actual  size  distribution  for  cases  2 and  3 haa  been  plotted  (Fig.  8}  according  to  Eq.  (20).  It  is  easily 
seen  that  over  Hamburg  the  smaller  particles  are  much  more  numerous  than  over  Edinburgh.  Between 
r§  ± 0.05  m*  31  percent  of  the  total  number  are  found  over  Hamburg,  but  only  15  percent  over  Edinburgh. 
It  is  safe  to  say  that  practically  all  particles  over  Hamburg  are  found  between  0.3  and  0.9  m «nd  between  0.5 
and  1.3  over  Edinburgh. 
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Fig.  7.  The  efficiency  factor  E\f  as  a function  of  the  size  parameter  for  m * 4/3  in  a logarithmic  scale  (after  Foitaik)  for 
case  3.  Curve  a for  h m 7.07;  curve  b for  h * 5.20;  curve  c for  h =-  3.87.  Curve  d;  measured  extinction  coefficient  (in 
relative  scale)  as  a function  of  wave  number  n (after  Wilson).  The  data  between  n * 2.4  and  2.6  are  considered  less  reliable 
by  the  author  (Wilson).  Curve  d (between  n = 1.6  and  2.4)  matches  curve  b best,  but  does  not  match  my  of  the  others. 


It  is  very  hard  to  draw  further  conclusion,  because  both  observations  are  based  more  or  less  on  one 
exposure  taken  at  a distinct  time.  If  a series  of  exposures  had  been  taken,  perhaps  a change  with  time  of  the 
size  distribution  may  have  been  observed. 

If  m « 4/3  is  not  the  appropriate  value  for  the  refractive  index,  Table  4 (b)  shows  the  numerical  results 
according  to  Table  3 (Section  4.13).  Under  the  headings  "Wellmann”  and  "Wilson,”  the  values  given  in 
their  publication  have  been  entered.  Under  the  heading  "Penndorf,”  those  calculated  according  to  Table  3 
are  given.  For  ice  particles  at  very  low  temperature  (m  = 1.25),  r0  is  much  larger,  namely,  0.660  m over 
Hamburg  and  1.08  fj,  over  Edinburgh,  but  the  total  numbers  are  smaller.  If  m = 1.46,  namely  for  wood  oils, 
rQ  = 0.380  p over  Hamburg  and  0.59  /x  over  Edinburgh  and  the  total  number  is  about  10*/ cm2  column.  Higher 
values  of  m are  unlikely  but  the  numerical  values  are  given.  The  effect  of  anisotropy  of  the  particles  will 
alter  ro  only  by  a few  percent  and  is  insignificant  in  comparison  to  other  uncertainties. 
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Fig.  8*  Size  distribution  for 
cases  2 and  3 and  m = 4/3. 
The  ordinate  gives  the  number 
of  particles  N contained  within 
a radius  range  of  0.1  The 
shaded  area  shows  the  number 
of  particles  contained  in  the 
range  rQ  ± 0.05  fi. 


In  general,  it  cap  be  stated  that  at  that  particular  time  the  radius  r0  was  between  0.38  m and  0.66  p over 
Hamburg  with  a total  number  Nt  between  3.7  X 107  and  1.26  X 10®/cm2  column,  whereas  over  Edinburgh 
the  numbers  for  r0  are  0.59  m and  1.08  m and  Nt  between  2.68  and  9.12  X 107/cm2  column.  Furthermore, 
there  existed  a spread  over  a size  range  which  changes  with  respect  to  the  m selected,  but  the  character  of 
Fig.  8 for  m * 4/3  will  remain. 


5.  ANOMALOUS  EXTINCTION  BETWEEN  3500  A AND  7000  A FOR  SEPTEMBER  1950 

Now  the  numerical  values  of  the  extinction  on  27  September  1950  can  be  calculated.  The  extinction  of 
the  normal  air  & is  taken  from  Linke  (1942,  p.  249)  for  a clear  day.  It  contains  the  extinction  due  to  molecu- 
lar scattering,  scattering  by  large  particles  which  contaminate  the  air  on  a "normal  clear"  day  and  absorption 


Table  5.  Extinction  coefficients  and  absolute  intensities  of  the  solar  radiation  for  the  blue  sun. 
Steps  100  A between  3500  A to  7000  A 
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Table  5.  Extinction  coefficient  and  absolute  intensities  of  the  solar  radiation  for  the  blue  sun  - — Concluded. 
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by  atmospheric  constituents.  Other  data  of  0 were  also  used;  however,  since  & contributes  only  a small 
amount  to  the  total  extinction,  such  differences  are  really  insignificant.  For  unit  air  mass  jS  decreases  stead- 
ily from  0.755  at  3500  A to  0.434  at  7000  A,  as  shown  in  Fig.  9 and  Table  5. 

The  extinction  of  the  smoke  particles  is  calculated  in  steps  of  100  A,  using  Eq.  (1)  and  the  thickness 
s = 3.7  km  for  cases  1,  2 and  3.  The  results  obtained  are  reproduced  in  Fig.  9 and  in  Table  5. 

The  anomalous  extinction  is  significantly  different  between  cases  1 and  2 on  one  side  and  case  3 on  the 
other.  Over  Hamburg  it  is  small  in  the  blue  and  red  regions,  and  reaches  a maximum  at  4900  A in  case  1 
and  6200  A in  case  2.  The  differences  for  ys  in  cases  1 and  2 arise  from  the  different  E functions  used; 
case  2 is  the  more  realistic  one.  Over  Edinburgh  a distinct  minimum  occurs  at  4300  A,  the  extinction  being 
large  in  the  ultraviolet  and  red  parts  of  the  spectrum.  In  all  three  cases,  ys  ^ 0:  for  3500  A it  is  greater  by  a 
factor  of  1,5,  2 and  5 for  cases  1,  2 and  3,  respectively,  and  for  7000  A by  a factor  of  3,  4 and  10  for  cases  1,  2 
and  3,  respectively.  Thus,  the  various  values  of  0 to  choose  from  did  not  alter  the  result.  Moreover,  slight 
changes  of  the  normal  atmospheric  transparency  between  the  two  days  Wellmann  and  Wilson  selected  are 
without  any  significance  to  the  result. 

The  narrow  region  measured  by  Wellmann  is  clearly  conceived  from  the  diagrams.  The  linearity  of  ys 
is  more  or  less  restricted  to  a small  interval  he  measured.  The  total  extinction  0 + 7*  is  also  plotted  in  Fig.  9 
and  the  numerical  results  are  given  in  Table  5.  For  Wellmann’s  observation,  both  cases  show  that  the  total 
extinction  is  small  in  the  blue  part  of  the  spectrum  and  the  maxima  occur  at  5700  A.  For  Wilson’s  observa- 
tion, the  total  extinction  shows  a broad  minimum  between  4200  A and  5600  A. 


6.  SPECTRUM  OF  COLORED  SUN 

6.1.  SPECTRUM  AND  INTENSITY  OF  A "BLUE”  SUN 

The  next  step  is  to  calculate  the  solar  radiation  on  the  ground  using  Eqs.  (2)  and  (3).  The  values  for  J0X 
are  taken  from  Linke  (1943);  the  extinction  coefficients  0 + ys  are  taken  from  Table  5 as  described  in 
Section  5,  and  thus  J^x r and  Ia\{  are  calculated  in  steps  of  100  A for  the  region  3500  A to  7000  A,  and  for  the 
zenith  distances  f — 53°,  61°,  67°  and  71°.  The  spectral  region  was  selected  for  two  reasons,  first,  the 
absorption  due  to  ozone  and  water  vapor  can  be  neglected;  second,  the  human  eye  is  sensitive  only  in  this 
region.  The  zenith  distances  selected  correspond  to  Dietze’s,  Wellmann’s  and  Wilson’s  observations. 

The  absolute  values  of  the  total  radiation  received  on  the  ground  are  very  small,  only  a fraction  of  what 
is  received  on  a normal  clear  day.  If  the  radiation  incident  on  the  atmosphere 

X- 7000  A 

Io  m I foxdX 

X- 3600  A 

and  1st a8  Wfll  a8  are  calculated  accordingly,  the  ratio  iAt/h  indicates  very  clearly  the  radiation  received 
on  the  ground.  The  numerical  results  as  found  in  Table  6,  show  that  the  solar  radiation  Jo  is  reduced  by  a 
factor  10“*  to  10“*.  Furthermore,  the  variation  caused  by  the  zenith  distances  is  extremely  large  (i.e.,  by  a 
factor  of  10  for  going  from  67°  to  53°).  Thus  it  is  seen  that  the  smoke  layer  reduces  the  solar  radiation 
tremendously,  comparable  to  a thin  water  cloud. 
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Table  6.  Direct  tolar  radiation  received  at  the  ground. 

The  ratio  Iff/h  corresponds  to  a normal  clear  day,  the  ratio  lj\/h  for  the  "blue  sun” 
observations,  and  the  ratio  1a*  / W for  Dietze’s  observation.  All  ration  refer  to  the  region 
3500  A to  7000  A only,  except  those  measured  by  Dietze. 


Zenith 

distance 

Computed 

Observed 

ls/U 

I All 0 

Ia'/W 

Case  1 

Case  2 

Case  3 

53° 

0.43 

0.03 

0.014 

5 6 X 10~< 

0.15 

61* 

0.35 

0.013 

0.0049 

9.4  X 10-6 

67° 

0.28 

4.7  X 10-< 

1.4  X 10-< 

1.0  X 10“* 

71° 

0.23 

.... 

1.2  X 10'« 

Table  7.  Ratio  of  relative  intensities  for  the  wavelength,  where  the  intensity  of  the  blue  sun  is  a maximum. 


Case 

1 

2 

3 

^max  iu  A 

4100 

4100 

4600  and  4700 

Zenith  distance 

U/h 

1a!  Ix 

Ja//o 

Ia/In 

Ia/U 

I a/ Is 

53° 

1.40 

1.69 

1.79 

2.16 

1.14 

1.21 

61° 

1.46 

.1.86 

2.06 

2.63 

1.18 

1.26 

67° 

1.65 

2.21 

2.46 

3.29 

1.22 

1.38 

71° 

.... 

.... 



1.30 

1.48 

Nore:  The  abbreviated  notation  I a/ h *•  used  for:  - I— —X  and  similar  for  I a/  f.v. 

I A ,6600/  f(),6600 


Dietze  observed  a larger  value  of  15  percent  for  Ia'/Iq  on  28  September.  The  reason  for  this  dis- 
crepancy of  our  values  and  the  observed  one  is  threefold.  First,  the  reduction  was  largest  on  27  September 
and,  therefore,  Nt  must  have  been  smaller  when  he  took  his  observations;  second,  the  actinometer  measured 
the  red  and  infrared  radiation  too,  which  is  excluded  from  our  calculations.  The  contribution  from  the  red 
part  of  the  spectrum  is  large.  Furthermore,  his  red  ratio  (Table  1)  indicates  this  fact,  which  is  due  to  the 
decrease  of  E with  wavelength  (part  (a)  of  Fig.  1).  Third,  the  actinometer  has  a large  angular  aperture  of 
several  degrees,  and  sky  radiation  around  the  sun  is  measured  also.  As  Dietze’s  data  show  (Table  1),  the 
sky  radiation  is  much  larger  than  normal  ( H/I  = 2.4  instead  of  0.09  on  a clear  day).  If  Ia/Iq  was  about  0.1 
in  Dietze’s  case,  Nt  may  be  assumed  to  be  3.5  X 107  (1/cm2  column).  Such  a variation  seems  very  reason- 
able and  the  total  amount  certainly  varied  (see  Section  6.3). 
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The  most  interesting  features  of  the  gpectrum  are  obtained  from  Fig.  10.  Here  relative  values  are 
plotted  using  Jo  **  I a In  = 1.0  for  X = 5500  A,  corresponding  to  the  maximum  sensitivity  of  the  eye. 

The  most  impressive  spectrum  is  certainly  that  of  Fig.  10  (b)  for  case  2,  but  the  outstanding  characteristic, 
the  maximum  in  the  blue,  is  the  same  in  all  three  cases.  The  intensity  in  the  blue  increases  with  increasing 
zenith  distance,  if  relative  values  as  defined  above  for  Fig.  10  are  used.  For  X ~ 4100  A,  the  relative  inten- 
sity is  1.8  times  Jox  for  f = 53°,  and  2.4  for  f = 67°  for  case  2 and  smaller  for  case  3 as  seen  from  Table  7. 

The  increase  in  the  yellow  and  red  parts  of  the  spectrum  starts  in  case  1 which  is  already  at  5700  A,  and 
in  case  2 at  6700  A.  Thus  case  2 is  certainly  more  realistic  from  this  point  of  view  also.  This  increase  is  a 
consequence  of  the  E function,  which  decreases  rapidly  from  the  first  maximum  towards  0 for  a = 0 or 
X — *oo  as  can  be  seen  from  Fig.  1,  as  well  as  from  Figs.  5 and  6.  Case  3,  Fig.  10  (c),  shows  a different 
behavior  because  the  maximum  in  the  blue  stands  out  much  clearer,  even  if  this  increase  is  smaller,  as  in 
cases  1 and  2 (Table  7).  However,  the  relative  intensity  in  the  red  part  of  the  spectrum  is  much  lower  than 
the  ratio  J0tX/  Jo,moo  ; at  7000  A I a*/ I a mm  » less  than  one-half  of  J0px/  Jo.wooi  where  as  in  case  2 it  is  nearly  equal. 

The  results  of  the  calculations  definitely  confirm  that  the  energy  distribution  of  the  solar  spectrum  has 
been  complete!)  altered  and  that  a maximum  exists  in  the  blue  region,  provided  that  the  radius  of  the  particles 
is  about  0.4  to  0.6  m and  that  the  total  number  is  very  large,  namely,  « lOVcm.  These  are  the  two  important 
conditions  which  cause  a "blue”  sun. 

In  the  various  general  descriptions  of  the  "blue”  sun  of  September  1950,  reports  are  also  given  of  the 
same  phenomenon  after  dust  storms  near  desert  regions.  In  such  cases  the  particles  are  probably  quartz 
dust.  The  refractive  index  of  quartz  is  m = 1.55.  This  value  is  not  much  larger  than  the  values  used  in 
this  calculation,  and  from  Table  4 it  can  be  concluded  that  the  particles  must  be  smaller,  about  0.3  to  0.5 
and  the  number  108/cm2  column.  If  these  conditions  are  fulfilled,  then  the  solar  spectrum  is  similar  to 
that  of  Fig.  10  and  the  sun  may  look  blue  in  color. 


6.2.  HYPOTHETICAL  SPECTRUM  OF  A "GREEN”  SUN 

A green  sun  or  moon  has  been  reported  very  rarely  in  the  literature.  One  British  observer,  however, 
states  that  on  26  September  1950  the  moon  was  rising  as  a normal  orange-yellow,  but  later  the  color  was 
something  between  a sea-green  and  a turquoise -blue. 

It  can  be  concluded  from  Fig.  2 that  the  radius  of  the  particles  must  be  close  to  1 p for  a minimum  of 
extinction  in  the  green.  Thus  a calculation  similar  to  that  of  the  blue  sun  described  above  was  carried  out, 
using  r0  = 0.963  ^ Nt  = 2 X 107/cm2  and  three  size  distributions.  Case  4,  r = ro,  uniform  size;  case  5, 
same  size  spectrum  as  in  case  2,  with  h = 7.07;  case  6,  size  spectrum  with  h « 10,  i.e.,  a much  narrower 
distribution.  The  zenith  distances  chosen  are  53°,  61°,  67°  and  75°.  The  results  of  the  calculations  are 
given  in  Fig.  11  for  two  zenith  distances.  The  general  trend  can  be  obtained  from  the  selected  results  as 
published  in  Table  8. 

A relatively  sharp  maximum  can  be  obtained  only  for  a very  narrow  distribution  of  particles.  The 
reason  for  this  is  that  the  broader  the  size  spectrum,  the  less  pronounced  is  the  minimum  in  the  E function 
(see  Fig.  1).  There  is  practically  no  minimum  if  h < 7 or  <x  > 0.1,  which  means  that  the  anomalous  extinc- 
tion is  more  or  less  independent  of  wavelength.  Such  a neutral  extinction  shows  no  color  effects.  The 
intensity  maximum  becomes  more  pronounced  with  increasing  zenith  distance.  Hence,  only  for  a sun  and 
moon  close  to  the  horizon  can  the  green  color  effect  be  observed.  With  the  sun  or  moon  high  above  the 
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Fig.  1 1 . Spectral  intensity  distribution  of  tbe  solar  radiation  between  3500  A and  7000  A,  in  relative  units  for  a "green  tun.” 

(a)  Refers  to  oases  lf  2 and  3 and  aetulb  distance  of  tbe  sun  f ■ 53*. 

(b)  Refers  to  cases  1,  2 and  3 and  aeoith  distance  of  tbe  sun  f ■ 67®. 

The  relative  intensity  1.0  is  ebosen  for  X — 5500  A. 

horizon,  the  color  must  return  to  the  normal  color  of  the  objects.  The  conditions  for  a green  sun  are  very 
severe: 

(1)  r **  0.9  to  1.0  m and  a very  narrow  site  spectrum. 

(2)  Nt  107 /cm*  column. 
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Table  8.  Relative  intensities  Iax^/Ia f«oo  for  a green  sun. 

Three  wavelengths  are  selected*  namely,  3500  A,  7000  A and  Xmax  • 

Constants  chosen  are:  ro  — 0.963  Nr  ■ 2 X 107/cm*  column,  {*  — 53°,  61°,  67°  and  75°. 


Case  4:  uniform  size 


A In  A 

53° 

61° 

f 

67° 

75° 

3500 

0.11 

0.075 

0.049 

0.016 

5300 

I.  08 

1.10 

1.11 

1.15 

7000 

0.30 

0.23 

0.17 

0.079 

Case  5:  size  spectrum,  h * 7.07 


X in  A 

53° 

61° 

r 

67° 

75° 

3500 

0.26 

0.22 

0.20 

0.13 

4900 

1.05 

1.02 

.... 

5100 

.... 

.... 

1.02 

5300 

— 

1.00 

7000 

0.73 

0.70 

0.69 

0.62 

Case  6:  size  spectrum,  h = 10 


XinA 

53° 

61° 

r 

67° 

75° 

3500 

0.23 

0.19 

0.15 

0.09 

5400 

1.01 

1.01 

1.02 

1.02 

7000 

0.56 

0.51 

0.46 

0.34 

For  quartz:  m ■=  1.55,  the  same  results  are  obtained  as  calculated  above,  if: 

(1)  r « 0.58  m 

(2)  Nt  ^ 5.6  X lOVcm2  column. 


6.3.  PHYSIOLOGICAL  IMPRESSION  OF  SUCH  SPECTRA 

What  subjective  impression  does  such  a spectral  distribution  of  solar  light  as  received  by  the  eye  produce 
in  the  brain  of  the  human  observer?  This  is  a physiological  and  no  longer  a physical  problem.  We  shall  deal 
only  with  an  observer  whose  color  sensitivity  is  normal. 

The  mechanism  of  color  vision  is  but  partially  understood  at  the  present  time.  The  most  important 
problem  with  regard  to  atmospheric  optics  is  the  interpretation  of  spectrophotometric  data  by  the  eye.  The 
theory  developed  by  Helmholtz  is  based  upon  the  fact  that  any  color  stimulus  can  be  matched  visually  by  a 
mixture  of  the  proper  amounts  of  three  arbitrarily  chosen  stimuli,  which  are  called  primaries. 


y (GREEN ) 
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Fig.  12.  Chrommcity  diagram 
for  a green  and  a blue  aun.  A 
description  of  this  diagram  ia 
given  in  the  text.  The  numbers 
refer  to  zenith  distances  of 
the  sun. 


0.1  0.2  0.3  0.4  0.5  0 6 0.7 

> X (RED) 

The  excitation  functions  of  the  eye  for  three  such  primaries  are  defined  by  the  International  Commission 
on  Illumination  for  an  equal-energy  spectrum.  But  the  relative  magnitudes  of  the  three  functions  have  been 
chosen  so  that  the  areas  under  the  curves  are  equal  for  mean  noon  sunlight.  The  graphical  representation 
of  color  is  best  carried  out  in  a color  triangle  using  the  three  excitation  functions,  or  in  an  ordinary  Cartesian 
coordinate  system,  as  shown  in  Fig.  12.  If  the  values  are  normalized,  with  x\  *=  red,  y\  = green 
and  z\  * blue,  according  to 

*1  yi  *1 

X\  + yx  + Zi  xt  + yt  + %\  *i  + yx  + *i 

it  suffices  to  plot  x and  y,  because  only  two  of  the  three  quantities  are  independent.  Each  point  in  this 
diagram  represents  a certain  color. 

On  this  diagram  point  C represents  the-  center  of  gravity  and  corresponds  to  the  solar  spectrum  as  seen 
on  the  earth.  Furthermore,  it  shows  the  locus  of  spectrum  stimuli  as  a series  of  circles  starting  in  the  lower 
left  corner.  These  are  the  points  for  pure  monochromatic  colors  between  4000  A and  6000  A.  A straight 
line  from  C to  any  of  the  spectrum  loci  includes  all  points  representing  the  stimuli  obtained  by  additive 


/ 

4lO0 

I? 


, 5300 


GREEN 


\ 


r CASE  I X 
"BLUE"SUN  < CASE  2 + 

^ CASE  3 o 

"GREEN”  SUN  • 

EQUAL  ENERGY ~ Be 
BLACK  BODY  6000*  K 


\ 5500 
\ 


\ 5600 

°s 

\ 


\ 5 700 
O 

\ 

\ 

\o58° 

YELLOW  \ 


\ 


\ 

\ 


^4050 

\ 

^ 4800 


\ 4750 
O 

\4700 

°\ 

\ 4500 
<V4000 
Q 1 


BLUE 


PURPLE 


\ 6000 
O 


**•  - * *'-.-*.%*+**'■  . I**?**  ;***■*„, 


39 


Table  9.  Trichromatic  coefficients  x%  y anu  z,  dominant  wavelength  (X^)  and 
colorimetric  purity  (p)  for  the  blue  and  green  solar  spectra  as  functions  of  aenith 
distance  f and  particle  size  spectrum. 

(a)  Blue  Sun  Case  1 


t 

53° 

61° 

67° 

* (red) 

0.327 

0.330 

0.333 

V (green) 

0.314 

0.314 

0.309 

z (blue) 

0.339 

0.356 

0.356 

Xrf 

4850 

P 

0.07 

negative 

negative 

color 

blue 

purple 

purple 

Case  2 

r 

53° 

61° 

67° 

x (red) 

0.286 

0.278 

0.261 

y (green) 

0.283 

0.273 

0.262 

z (blue) 

0.430 

0.449 

0.477 

U 

4750 

4700 

4750 

P 

0.21 

0.24 

0.32 

color 

blue 

blue 

blue 

Case  3 

r 

53° 

61° 

67° 

71° 

* (rod) 

0.282 

0.274 

0 265 

0.257 

y (green) 

0.313 

0.309 

0.305 

0.299 

z (blue) 

0.405 

0.417 

0.429 

0.444 

\d 

4850 

4850 

4850 

4850 

P 

0.18 

0.22 

0.25 

0.28 

color 

blue 

blue 

blue 

blue 

(b)  Green  Sun 

zenith  distance  * 67" 

Case  4 

Case  5 

Case  6 

x (red) 

0.308 

0.336 

0.333 

y (green) 

0.415 

0.355 

0.389 

s (blue) 

0.277 

0.309 

0.278 

x< 

5350 

5600 

5550 

P 

0.17 

0.05 

0.16 

color 

green 

yellowish-green 

yellowish-green 
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mixture  in  various  portions  of  the  extremes  of  the  spectrum.  Therefore,  it  is  called  a mixture  (or  chro- 
macity)  diagram. 

There  are  other  factors  besides  chromacity  which  determine  the  color,  e.g.,  hue,  saturation  and  brilliance. 
Three  objective  parameters  are  chosen  to  correspond  in  a general  way  with  these  physiological  attributes. 
They  are: 

(1)  dominant  wavelength 

(2)  colorimetric  purity 

(3)  brightness  (or  luminosity). 

The  dominant  wavelength  is  defined  as  the  wavelength  where  a straight  line  between  C and  the  color 
point  intersects  the  curve  for  pure  monochromatic  radiation  (spectrum  locus).  The  colorimetric  purity  is 
the  ratio  of  the  luminosity  of  the  monochromatic  radiation  to  the  total  luminosity.  The  luminosity  is  simply 
the  sum  of  the  luminosities  of  the  monochromatic  radiation  and  white  radiation. 

These  considerations  are  now  applied  to  the  spectra  calculated  in  Sections  6.1  and  6.2.  The  numerical 
constants  and  formula  are  taken  from  Judd  (1933)  and  the  results  are  condensed  in  Table  9 and  plotted  on 
Fig.  12.  In  the  case  of  the  "blue”  sun,  case  1 gives  values  for  the  uniform  particle  size  which  are  very  close 
to  the  solar  color  (crosses  in  Fig.  12).  The  purity  is  0.07  for  f = 53°;  for  the  other  two  the  color  should 
be  white  with  a purplish  hue.  Thus  this  poor  result  shows  that  Wellmann’s  data  cannot  be  interpreted  by 
particles  of  uniform  size.  But  in  case  2 the  values  are  better;  the  dominant  wavelength  is  about  4750  A, 
which  is  blue.  The  colorimetric  purity  is  0.2  to  0.3.  Finally,  the  results  for  case  3 resemble  those  for 
case  2;  the  dominant  wavelength  corresponds  to  a blue  color  and  the  colorimetric  purity  is  0.18  to  0.28. 
The  colorimetric  purity  increases  with  increasing  zenith  distance.  The  color  should  not  be  pure  blue,  but  a 
rather  bluish  white,  because  the  colorimetric  purity  is  0.2  to  0.3  for  the  zenith  distances  in  question.  From 
these  data  it  is  seen  that  the  sun  or  moon  is  blue  near  the  horizon  and  becomes  more  and  more  the  normal 
yellow  - white  the  higher  it  is  in  the  sky.  It  is  also  seen  that  the  color  effect  only  occurs  for  particles  having  a 
size  spectrum  and  not  for  those  of  uniform  size. 

For  the  green  sun,  the  results  of  the  calculations  are  summarized  in  Table  9 (b)  and  in  the  mixture 
diagram.  The  values  for  zenith  instances  f * 67°  show  that  an  effect  only  occurs  for  very  uniform  particles; 
the  colorimetric  purity  is  0.17  for  uniform  particles,  for  a relatively  small  spread  of  h « 1C  it  is  0.16,  but  for 
a wide  distribution  with  h « 7 it  is  only  0.05.  The  color  of  the  pure  spectrum  locus  is  green  in  all  three 
cases.  The  colorimetric  purity  also  increases  with  increasing  zenith  distance,  as  for  the  blue  sun. 

The  final  computations  show  very  clearly  that  the  physiological  impression  of  the  spectra  calculated  in 
Section  6.1  is  a bluish-white.  Furthermore,  the  chromacity  diagram  proves  that  Wellmann’s  and  Wilson’s 
physical  measurements  confirm  the  numerous  observations  of  a blue  sun  during  September  1950.  It  has 
been  argued  by  Elvey  (1934)  that  the  appearance  of  a blue  sun,  during  the  observations  of  a duststorm  in 
1934  was  a "psychological”  and  not  a "physiological”  effect,  and  was  dependent  on  the  change  of  contrast 
between  sun  and  sky.  Since  he  used  only  a blue  and  yellow  filter  and  not  spectroscopic  measurements  — 
as  did  Wellmann  and  Wilson  — his  conclusions  are  based  on  insufficient  physical  data.  This  study  endeav- 
ors to  show  distinctly  that  very  precise  spectroscopic  observations  are  required  to  deduce  such  data  as  given 
in  Sections  4,  5 and  6. 

These  studies  confirm  that  the  observers  were  correct  in  describing  a "blue”  sun.  This  phenomenon 
was  not  fictitious,  but  real.  It  occurs  very  rarely,  simply  because  the  conditions  for  its  occurrence  are  very 
severe;  namely,  the  radius  must  be  around  0.5  g,  the  size  distribution  must  be  narrow  and  the  total  number 
larger  than  10*/cm*  column. 
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